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ABSTRACT 

We use statistical equilibrium equations to investigate the IRAC color space of shocked molecular 
hydrogen. The location of shocked H2 in [3.6] — [4.5] vs [4.5] — [5.8] color is determined by the gas 
temperature and density of neutral atomic hydrogen. We find that high excitation H2 emission falls 
in a unique location in the color-color diagram and can unambiguously be distinguished from stellar 
sources. In addition to searching for outflows, we show that the IRAC data can be used to map the 
thermal structure of the shocked gas. We analyze archival Spitzer data of Herbig-Haro object HH 54 
and create a temperature map, which is consistent with spectroscopically determined temperatures. 
Subject headings: ISM: jets and outflows — ISM: lines and bands — ISM: individual(HH 54) - 
molecular processes 



1. INTRODUCTION 

Protostellar outflows have strong molecular hydro- 
gen emission in the wavelength range covered by the 
Spitzer Infrared Array Camera (IRAC). Spitzer stud- 
ies of known outflows reveal that shocked H2 emission 
appears particularly strong in the 4.5 p,m IRA C band 
(jNoriega-Crespo et al.|[200l iTeixeira et al.ll2008l ). Many 
groups are beginning to visually inspect IRAC data to 
search for outflows and objects with extended H2 emis- 
sion based on the stro ng emission in the 4.5 /im band. 
iSmith fc Rosen! (|2005h created synthetic Spitzer images 
from their models of precessing protostellar jets. These 
models calculated the population of the first three vi- 
brational levels by solving for statistical equilibrium and 
assuming local thermal equilibrium (LTE) for the rota- 
tional levels. Their simulations showed that the emission 
in the 4.5 /im band to be the strongest, which was consis- 
tent with observations. In taking a ce nsus of the young 
stellar objects (YSOs) in NGC 1333, iGutermuth et all 
(2008) empirically determined an IRAC color cut based 
on observations of known shocked emission within NGC 
1333. This was used to remove any possible shocked 
emission in their source list of YSOs. Due to the multi- 
tude of lines in the IRAC bands it was thought that infor- 
mation on the physical parameters of the gas could not 
be ascertained from the Spit zer IRAC data. I n ana lyzing 
the IRAC images of IC 443. iNeufeld fc Yuan! (f2008h cal- 
culated the IRAC band ratios for shocked H2 using the 
13 strongest lines covered by IRAC but only included 
collisional excitation by H2 and He in their calculations. 
They found the measured flux in the 3.6 /im band to be 
stronger than predicted in their calculations, which may 
have been due to neglecting collisional excitation with 
atomic hydrogen. Until now, the color space of shocked 
H2 emission due collisional excitation with atomic hydro- 
gen, He, and H2 has not been calculated. 

In order to use Spitzer IRAC data to find outflows and 
study their properties we have calculated the IRAC color 
space of shocked H2 emission. 

2. CALCULATIONS 
Electronic address: jybarra@astro.ufl.edu 



We have calculated the intensities of shocked molecu- 
lar hydrogen and determined the location of the shocked 
emission within IRAC color space. At typical shock tem- 
peratures and densities, the excitation of molecular hy- 
drogen is through collisions with H atoms and He atoms, 
and with ground state ortho- and para-H2 molecules. 
The typical densities in outflows are less than critical 
so we do not assume LTE. Instead the populations of the 
first 47 ro- vibrational excited states were calculated from 
solving the equations of statistical equilibrium where 
we set n(He)/nH = 0.10, njj = n(H) + 2n(H2), and 
n(H)/n(H2) = 0.10. The atomic hydrogen fraction was 
set to the median value consistent w ith shock models 
and the temperature rang e we chose (jLe Bourlot et al.l 
l2002t iTimmermannl I1998D . To solve the equations we 
employe d a n on-LTE code based on the method by 
iLi et all (|1993f) . We used the latest H-H 2 non-reactive 
collisi onal rate coefficients calculated bvfWrathmal l et ail 
(|2007f h The reactive collis ional rate coefficients we re de- 
rived from the relations of iLe Bourlot et all (jl999). The 
rate coefficients for He- H? and H2-H2 collisions are from 
ILe Bourlot et all (jl999fh The de generacy of the states is 
given by gj — (2 J+ 1) for even J, and gj = 3(2 J+ 1) for 
odd J. The quadrupole trans ition probabilities used are 
from fWolnicwi cz et al.l (|1998l ) . We calculated the pop- 
ulations of the states over a wide range of atomic hy- 
drogen densities, n(H) = 50 — 10 5 cm~ 3 , and gas tem- 
peratures, T = 1000 — 6000 K. The relative intensities 
of 49 lines that fall within the IRAC bands (3.08 /jm 
< A < 10.5 /Ltm) were determined and used to calculate 
the IRAC band fl uxes using the la test published IRAC 
spectral response ( Hora et aLll2008| ) and calibration data 
l|Reach et al.l 12005 ). Table 1 lists the fractional contri- 
bution of the strongest H2 lines to each IRAC band for 
ra(H) = 10 4 cm" 3 and temperatures 2000 K and 4000 K. 

3. H 2 EMISSION IN IRAC COLOR SPACE 

Figure 1 shows location of shocked H2 gas in IRAC 
[3.6] - [4.5] vs. [4.5] - [5.8] color space. We find that the 
observed emission in the bands is a function of kinetic 
gas temperature and atomic hydrogen density. We re- 
strict our plot to gas temperatures below 4000 K where 
H2 emission is expected to be the dominant. Shocks that 
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TABLE 1 

Fractional contribution of the strongest H2 lines to the 
IRAC bands 8 - 



Transition A (jjsa) IRAC T=2000 K T=4000 K 

0.20 
0.08 
0.06 
0.02 
0.10 
0.13 
0.34 
0.04 
0.29 
0.12 
0.07 
0.33 
0.05 
0.14 
0.51 
0.10 
0.13 
0.66 
0.11 
0.12 



a fractional contribution to the total emission from the H2 lines in 
the bands convolved with the IRAC spectral response for n(H) = 
10 4 cm" 3 



1 i r 




-1012 

[4.5H5.8] 



Fig. 1.— IRAC [3.6] -[4.5] vs. [4.5] -[5.8] color-color plot indicat- 
ing the region occupied by shocked H2 . Constant density (dotted) 
and temperature (solid) lines are indicated. The region to the left 
and above the da shed-dotted line was empirically determined by 
Gutcrmuth ct al. (2008) to contain outflows. 

would produce gas temperatures in excess of 4000 K are 
likely to be dissociative (J-type) decreasing the abun- 
dance of H2 molecules. These shocks are also likely to 
produce vibrationally excited CO emission and fine struc- 
ture [Fe II] emission in addition to the H2 emission. The 
CO molecule, which has a higher dissociation energy than 
H2, is able to survive at higher shock velocities and tem- 
peratures. In these high energy shocks, CO becomes vi- 
brationally excited and emits in v = 1—0 (4.45 /im < A < 
4.9 /im) lines and can contribute significantly to the total 
emission from the shocked ga s. (Gonzal ez- Alfonso et al.l 
2002; Dr aine fe Robergdll984h . Furthermore, the major- 
ity of [Fe II] lines in the range covered by IRAC fall within 
the 4.5 /im band. Consequently, 4.5 /im emission in ex- 
cess to the color space defined by H2 at 4000 K is likely 
to trace gas with T > 4000 K placing shocked emission 
into the upper left portion of the color-color diagram. 
We chose not to use the 8 /im band as there may be 



dust and PAH contamination in this band. PAH emis- 
sion is particularly strong in the 8 /im band and it is 
still unclear what contribution PAHs may have in the 
emission of shocked gas from outflows. Dust continuum 
emission also becomes a possibility as dust may survive 
the shock. Furthermore our calculations are restricted to 
a single temperature along the line of sight. In the case 
of multiple temperature components, the cooler gas will 
significantly contribute to the lower excitation pure rota- 
tional lines, 0—0 <S'(4) and 0—0 5(5), that dominate the 
emission in the 8 /im band. The 8 /im band would then 
trace the cooler temperature component while the other 
bands would be more sensitive to the hotter gas. Thus 
our analysis can include gas containing multiple temper- 
ature components along the line of sight but will still be 
restricted to analyzing only the hotter gas. 

As seen in Figure 1 , the shocked H2 lies in a well defined 
location in the color-color diagram. For comparison, the 
location of YSOs is also shown in th e figu re based on the 
criteri a of iGutermuth et all (|2008l ) and Megeat h et alJ 
(2004). Shocked H 2 gas with sufficiently high tempera- 
ture and atomic hydrogen density is found in a unique 
location on this diagram and can be distinguished from 
YSOs. This is co nsistent with the empir ically deter- 
mined color cut of IGutermuth et al.l (|2008l ). Therefore 
these colors offer an unambiguous method for searching 
for shocked gas from outflows/jets. However, there is 
overlap in IRAC color space between low temperature 
shocked H2 gas and protostellar sources. Consequently 
surveys searching for outflows using IRAC color analy- 
sis will be restricted to finding flows containing higher 
excitation gas. Additional data at different wavelengths 
(2MASS, MIPS, etc.) and/or morphology may be able 
to break this degeneracy. Our results are consistent 
with empirical observations of strong 4.5 /im emission 
in outflows and with the hydrodynamic simulations of 
iSmith fc Rosenl (f2005h . 

Using this color analysis, two methods can be applied 
to identify outflows in the images: 1) Visual inspection 
of 3-color images constructed out of the 3.6 /mi, 4.5 /im, 
and 5.8 /im band data with appropriate scaling to en- 
hance the shocked emission, and 2) analysis of the pho- 
tometry of the field where features consistent with colors 
of shocked H2 emission are selected. This can be accom- 
plished by evaluating the color pixel by pixel across the 
field. 

The location of the shocked H2 in color space depends 
on its temperature. Therefore color analysis provides 
a new way to probe the temperature structure of the 
gas. This can be accomplished by evaluating the colors 
pixel by pixel across the field. The colors can then be 
compared to the colors of shocked H2 and the tempera- 
ture of the gas can be estimated. Resulting temperature 
maps are restricted to temperatures between 2000 — 4000 
K and high atomic hydrogen densities. The color space 
for H2 at temperature less than 2000 K moves further 
into the color space of Class 1/0 objects and it becomes 
more likely to misidentify scattered light from YSOs as 
H2 emission. Note that, at low atomic hydrogen den- 
sities the constant temperature lines start to converge 
and thus temperature estimates from this region of color 
space may have large uncertainties. 

4. EXAMPLE 
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Spitzer IRAC detection and analysis of shocked H2 
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We applied our analysis to Spitzer archival data of the 
known outflow HH 54. We evaluated the IRAC colors at 
each pixel in the field of HH 54. The median value of 
the image was used to estimate the background and was 
subtracted from the images. Figure 2 shows the color- 
color diagram for the k nots (A,B,C,E,K,M) o f shocked 
H2 previously studied bv lGiannini et al.l (|2006( ). We plot 
all the pixels on the diagram that are encompassed by 
the knots. The majority of points fall within our calcu- 
lated color space for shocked H2 emission. One excep- 
tion is knot A which contains several pixels with colors 
that fall more than 3tr outside the 4000 K boundary. 
This region of color space with excess 4.5 /zm emission 
is consistent with higher gas temperatures and possible 
a dditional emissi o n from CO v = 1 — and [Fe II] . 

IGiannini et al.l (|2006f ) obtained spectroscopic data of 
various knots in HH 54 in the near-infrared and used 
H2 emission lines to estimate the rotational and vibra- 
tional excitational temperatures of the gas. We created 
a temperature map of HH 54 (figure 3) by selecting the 
pixels whose colors fall within the range we identified 
as belonging to shocked H 2 (T = 2000 - 4000 K) and 
[4.5] — [5.8] < 1.5, and then estimated the temperature 
of the gas based on the pixels location in color space. 
The estimated temperatures are compared to the vi- 
brational (v > 1) ex citation temperatures obtained by 
IGiannini et al.l (|2006l ) in table 2. Because typical atomic 
hydrogen densities in shocks are less than critical, the 
gas cannot be assumed to be in LTE. In these environ- 
ments the rotational temperatures are often far below the 
kinetic gas temperature, while the ro-vibrational v > 1 
excitation temperatures are close to the kinetic gas tem- 
perature. For most of the knots the temperatures we 
estimate from color analysis and the spectroscopically 
determined temperatures are consistent. 

There is an additional knot, labeled I 
(0^2000 = 12:55:54.8, £j2ooo=-76:56:06), 5 arcseconds 
to the east of knot C seen in the IRAC images that 
is not found in the NIR images of HH 54 which we 
identify as the mid- infrared count erpart to the optical 
knot HH 541 (jSandell et alj 119871) . The temperature 
of the gas in this knot peaks at 3500 K as seen in 
figure 2 and 3. This knot is also found on the edge of 
the [N e II] 12.8 ^m map of HH 54 by iNeufeld et all 
(2006). The presence of the fine-structure [Ne II] line is 
consistent with the high temperature structure of this 
knot. Knot B (T ~ 2600 K)is also found to be spatially 
coincident with strong [Ne II] 12.8 /im e mission. The 
distri bution of the [Fe II] 26 /im line (|Neufeld et al.l 
2006) covers a region containing knots A, B, M. We 
find that spatial distribution of fine-structure emission 
lines is consistent with regions within our temperature 
map where T > 2600 K. Knots A, B, and I are spatially 
aligned with each other as indicated by the green line in 
figure 3. These high temperature knots may trace the 
jet component of the outflow. This line points in the 
direction of IRAS 12496-7650 which is believed to be 
the source of the outflow. 

5. CONCLUSIONS 

We have quantitatively shown that analysis of Spitzer 
data can be used to discover and characterize emission 



from protostellar outflows. Shocked H2 with sufficiently 
high temperature and neutral atomic hydrogen density 

TABLE 2 
Temperature estimates of HH 54 



Knot 


q(J2000) 


<5(J2000) 


n~> a 

± go.a 


rp b 


A 


12:55:49.5 


-76:56:30 


3300±500 K 


3000±150 K 


13 


12:55:51.1 


-76:56:21 


2600±300 K 


3000±140 K 


C 


12:55:53.1 


-76:56:06 


2200±100 K 


2960±150 K 


E 


12:55:53.8 


-76:56:23 


2000±100 K 


2500±90 K 


K 


12:55:54.3 


-76:56:25 


2100±100 K 


2500±90 K 


M 


12:55:51.6 


-76:56:13 


2800±400 K 


3000±140 K 



a The estimated temperature is the average temperature corre- 
sponding to the pixels within each knot based on its location in 
IRAC color-color space. The uncertainty is the standard devia- 
tion of the temperatures corresponding to the individual pixels. 
b Vibrational excitation temperatures from Giannini et al. (2006) 
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Fig. 2. — IRAC color-color plot for identified knots in HH 54. 
Each point represents the colors at an individual pixel within the 
knots. Different knots are represented with different symbols. 



can be distinguished unambiguously from stellar objects 
in IRAC color space. IRAC color analysis is useful for 
studying intermediate-excitation shocked gas within the 
temperature range T = 2000 — 4000 K. Higher tempera- 
ture gas may contain significant contribution from ionic 
lines like [Fe II] and also from the CO v = 1 — emis- 
sion band. Even with this limitation, Spitzer IRAC data 
provides a useful tool in the study of outflows. In par- 
ticular, IRAC color analysis can be used to probe the 
thermal structure of the gas without the need of using 
spectroscopic data. 
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Fig. 3.— a) IRAC 4.5^m imag e of HH 54. Image ce nter is at a(J2000) = 12 h 55 m 49?5 <5(J2000) = -76°56'22". Individual knots are 
identified and labeled according to Giannini ct al. (2006). Knot I, not seen on the NIR, is identified as the optical counterpart to HH 541 
b) Temperature map of HH 54 based on IRAC color-color analysis for 2000 K < T < 4000 K. The black contour levels are 2000 K, 2500 
K, 3000 K, and 3500 K. The green line connecting the higher temperature knots A, B, and I points toward the proposed source IRAS 
12496-7650. 
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